Difficulties with the QPOs Resonance Model 

Paola Rebusco 

Kavli Institute for Astrophysics and Space Research, MIT, Cambridge, MA 02139, USA 

"To know that we know what we know, 
and to know that we do not know what we do not know, 

that is true knowledge." 

Nicolaus Copernicus 



Abstract 

High frequency quasi-periodic oscillations (HFQPOs) have been detected in microquasars and neutron star systems. 
The resonance model suggested by Kluzniak & Abramowicz (2000,2001) explains twin QPOs as two weakly coupled 
nonlinear resonant epicyclic modes in the accretion disk. Although this model successfully explains many features 
of the observed QPOs, it still faces difficulties and shortcomings. Here we summarize the aspects of the theory that 
remain a puzzle and we briefly discuss likely developments. 
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1. Introduction 

Many galactic black hole and neutron star sources 
in low-mass X-ray binaries show quasi-periodic 
variability in their observed X-ray fluxes. Some of 
the quasi-periodic oscillations (QPOs) are of the 
order of few hundred Hz (hence the name high fre- 
quency) and often come in pairs of twin peaks in the 
Fourier power spectra (see Ivan der Klis fc Murdinl . 



2OOOI : iRemillard fc McClintockll2006l . for a review) 
High frequency QPOs draw much interest be- 
cause their frequencies lie in the range of orbital 
frequencies few Schwarzschild radii outside the 
central source. Furthermore for different sources 
they scale with 1/M, wh ere M is the mass of the 
central compact object ( McClintock fc Remillardl . 
En06 ). These two characteristics make QPOs at- 
tractive tools to possibly test General Relativity 
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in strong fields. Many models have been pro- 
posed in order to understand HF QPOs: some of 



them involve orbital motions 



Kluzniak et al 



'1990'; 'Stella & Vietri', Iggs"; 'La mb fc Miller] . 12003 



accretion disk oscillations (e.f 


I.. Waffoner et al. 


2001': 'Kato*, '2001'; 'Rezzofla et al. 


, 2003: Katol. 2007 


Tassev & Bcrtschinger, 20071). None is definitive 



Kluzniak fc Abramowicz (2000,2001) proposed that 
certain properties of HFQPOs in neutron stars could 
be explained by resonant motions of accreting fluid 
in strong gravity, and predicted that in black hole 
systems two HFQPOs in a frequency ratio of 1:2, 1:3 



n:m, could be detected. Abramowicz fc Kluzniak 



( 200ll ) observed that t he newly discovere d HFQ- 
POs in GRO J 1655-05 ijStrohmaveil l200ll ) were in 
a 3 : 2 ratio, and proposed a more specific model 
of res onance. This idea g uided s uccessive s tudies 
(e.g., Abra mowicz et al.. 2003b : Rebuscol 2004; 
2004 



Horak et al 



Kl uzniak et al.l . 120041 ) aimed 



at clarifying and expanding the theory. Wlodek 
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Kluzniak talked about the accomplishments of the 
model (see the contribution in the same Proceed- 
ings), here we focus on its limits. 

2. High Frequency QPOs and NonUnear 
Resonances 

The fundamental features of the nonlinear reso- 
nance model are that: 

- HFQPOs arise from a nonlinear resonance in ac- 
cretion disks in General Relativity. 

- The frequencies of nonlinear oscillations depend 
on the amplitudes of the same oscillations. As a 
consequence the observed frequencies {v*) may 
differ from the eigenfrequencies (i^o) Sind can 
vary in time {v* =1^0 + Ai/{t)). We refer to this 
well-known nonlinear phenomenon as frequency 
correction. 

- The two frequencies of resonant modes are ap- 
proximately in the ratio of small integers, most 
likely in the ratio 3:2. 

3. Difficulties with the Nonlinear Resonance 
Model 

All models of QPOs are essentially dynami- 
cal models, that miss any emission mechanisms 
and any connections to the spectral states of the 
sources. The only semi nal steps in this direc- 



tion have been done by iBursa et al.l (|2004[ ) and 



Schnittman fc BertschingeT i 2004 ) ."who considered 



general relativistic ray-tracing from an oscillating 
slim torus and an orbiting hot-spot respectively. In 
what follows we concentrate on the difficulties in 
the pure dynamics, since the radiation problem still 
has to be addressed by all the models. 
HFQPOs from black holes are stable (for many 
years), while neutron star QPOs are variable and 
present a rich phenomenology: consequently the 
two classes of systems will be treated separately. 

3.1. Difficulties in Black Hole Systems 

Four microquasars display pairs of HFQPOs, 
whose centroids ratio is consi stent with 3 : 2 
([Remillard fc McClintoc^ . l2006l and see Figure [T|). 
Sometimes the twin peaks are observed simultane- 
ously, most often only one is present. GRS 1915+150 



has a second pair dStrohmaveil 2001 ) i n a 5 :3 fre- 
quency ratio ( Kluzniak fc Abramowica . 120021 ) . 
Contrarily to what it may seem at first sight, all 
these features are naturally explained in the res- 
onance model. Indeed the absence of a peak may 
be due to the relativistic modulation of the emit- 



ted light (jBursa et al.l . 120041 ). Moreover although 



the strongest resonance is expected to occur when 
the frequencies are in 3 : 2 ratio, lower order reso- 
nances (such as 5 : 3) and sub harmonics can still 



be excited and become visible ( Abramowicz et al 



2003bl : iRebuscol . |2004 iHorak et all 12004 ). In the 



resonance picture the vicinity of the observed ratio 
to commensurate ratio indicates that the frequency 
corrections are small. 

The real problems of the model are more substan- 
tial: how are the QPOs excited? How are they 
couple d? In a preliminary toy-model b ased on the 
idea of lAbramowicz fc KluzniakI (|200l[ ). the motion 
of a single particle on a perturbed geodesic was 
analyzed. In t his simple model an ad-hoc parame- 
ter w as added ( Abramowicz et al. ■ l2003bl : [RebuscoL 
|2004) . that would provide the energy to increase 
the amplitudes of the resonant oscillations. This ap- 
proach works well mathematically, but the physics 
remains unexplained. A more detailed fluid model 
of the disk should be developed in order to jus- 
tify this excitation mechanism and to understand 
whether for example a form of disk instability (e.g., 
IPapaloizou fc Pringlel 1984 . or magnetic instabil- 
ity) could arouse the oscillations. 



Accretion disks are highly turbulent. IVio et al 
([2OO6) investigated the possibility that stochastic 
turbulence could excite and feed the resonance. To 
this purpose a white noise term was added to the 
vertical component of the perturbed geodesic and 
it was found that a small noise could indeed trigger 
the epicyclic resonance (see figure [2|) . This result 
opens an interesting perspective. Still a more real- 
istic scenario, that takes into account colored noise 
and stochastic forcing of a disk (rather than a single 
particle) should be evaluated. 

Finally in the alpha of the QPOs resonance theory 
it was always mentioned that gravity and pressure 
would provide the coupling between the resonant 
modes. However it was recently estimated that in 
thin disks and slender tori such a coupling is too 
weak to actually transfer energy between epicyclic 
modes in a reasonably short time (Jifi Horak, 
private communication). This difficulty holds in 
neutron star systems as well. One possibility could 
be to consider the coupling of epicyclic modes in 
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Fig. 1. Lower (fdomn) versus upper (uupp) frequencies in observed twin peak HFQPOs. The four microquasars lie on the lower 
left part of the plot, on the line with slope 3:2. Neutron star systems (all the others) cross the 3 : 2 line, but at different 
times of observation may be in a different position and ratio (courtesy of Michal Bursa). 



different geometric configurations. Still this may 
originate new dilemmas: for instance, in a thick 
torus the epicyclic oscillations (if not damped) do 
not simply scale with the mass of the central object, 
but depend st rongly on therrnqdyna mical quanti- 
ties. Moreover Sramkova et al.l ( 2007t ) showed that 
in a symmetric newtonian slender/thick torus there 
is no pressure coupling between epicyclic modes. 
Non-axisymmetry may play an important role. 
Another option would be to consider non-epicyclic 
modes, but one should understand why two non- 
epicyclic modes are more important than others. 
Although much has been done, the excitation and 
couphng of modes still needs to be developed and 
a close interaction between analytic and numerical 
work is fundamental. Simulated axisymmetric tori 
slightly off equilibrium do not incontrove rtibly pro- 
duce QPOs : they need to be "kicked" (iLee et al" " 



20041 : IZanotti et all 120051 : iBlaes et all . |2006(). Once 



the modes are excited (via velocity, density per- 
turbations or periodic forcing), then a spectrum of 



modes arise, with some modes in 3 : 2 ratio. One 
of the challenges consists to identify the simulated 
modes with analytic solutions and to tag if there is 
any resonances. 

We conclude this section on rn i croqu asars pointing 
at another puzzle: iTorok et al. I (|200d) reported that 
the spin estimates from the resonance model do 
not match the estimates obtained by fitting spec- 
tral continua. This is a puzzle that any relativistic 
QPOs model must address. 

3.2. Difficulties in Neutron Star Systems 

As we have seen there is no general agreement 
about the nature of HFQPOs in black hole systems. 
For neutron stars the situation is even more complex 
and the scientific community is torn between those 
who think that HFQPOs in neutron star systems 
have the same origin as in microquasars, and those 
who sustain that it is a different phenomenon. In 
the scenario proposed by EE ramowicz fc KluzniakI 
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Fig. 2. Spectra (upper panels) and phase diagrams (lower panels) for radial (p) and vertical (z) perturbed geodesies. The 
displacements are in units of Schwarzschild radii, the frequencies are scaled to kHz (assuming a central mass M of 2 Mq). 
The left plots are noise free, while the right have a small noise in the z direction. Clearly the amplitudes of oscillation are 
much greater in presence of stochastic forcing. Note that if turbule nce becom es too strong, then any periodicity is lost and 
the motion becomes chaotic. See section [3. II and the original paper l lVio et al.L 2006) for a more extensive explanation. 



([20011) HFQPOs have the same root in both classes 
of sources, with some important caveats. 
Twin HFQPOs are observed in accreting neutron 
stars: however both peaks can be in a different po - 
sition at a different time (e.g.. Ivan der "Kii3 . l2005l) . 
The frequency shift can be higher than a factor of 2 
(see Figure[T]). The ratio of the two centroids is some- 
times in the vicinity of 3 : 2, other ti mes it differs 
from the exact ratio (up to ~ 10%). iBelloni et al 



(j2005f ) argue that in neutron s tar sources there is 
not a preferred fixed ratio (as lAbramowicz et aT 



2003al . claim) and the occurrence of ratios close to 
3 : 2 is biased. 

In the resonance model these QPOs are still inter- 
preted as resonant epicyclic modes, with amplitudes 
larger than in the case of black holes. The neutron 
star surface plays an important role: the X-ray emis- 
sion is indeed modulated at the boundar y layer , 
where the accreted matter hits the star ( HorakL 



2005t lAbramowicz et al.l . 120071 ). The divergence 
from the exact ratio and the possibility of more than 
one ratios (as suggested bv lBelloni etall . l2005f ) are 
natural consequ ences of the nonlinear nat ure of the 
oscillators (e.g.. Landau &: Lifshitz . 1969t ) 
The bright Sco X-1 is often taken as a paradigm 
for neutron star sources. Different detections lie 
on a line with a slope close , but not equal to 



3 : 2 dyan der Klis et al 



19971) . lAbramowicz et al 



(|2003bl ) and iRebuscol (|2004) reproduced part of 



this line, by considering the frequencies and ra- 
tio shifts as frequencies corrections (Ai/(i)) due 
to nonlinear resonance. This result is qualitatively 
interesting, but the toy-model they adopt is too 
simple. Moreover in the nonlinear theory the lowest 
order frequency corrections scale with the square of 
the amplitude of the perturbation. Shifts as large 
as those observed would mean so high amplitudes 
that the adopted expansion would fail. 
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HFQPOs in neutron star sources may be related 
to the spin of the central object. In the very first 
detections the difference of the frequencies of the 
two peaks was found to be consistent with being 
constant and equal to the spin of the neutron star. 
When more sources were observed, they were di- 
vided into two classes: slow {vspm < 400 Hz) and 
fast rotators, in which the difference was equal to th e 



spin or half of it respecti velv (,Muno et al.l. 120011). 



From the study of Sco X-l lAbramowicz fc KluznialJ 
(|2004h pointed that HFQPOs have "little to do 
with the rotation of a stellar surface or any mag- 
netic fie ld structure ancho r ed in the star". Very 
recently iMendez fc Bellonil (|2007t ) re-examined all 
the available data and found that the difference of 
HFQPOs is approximately constant but not related 
to the spin and this issue is currently under debate 



(jStrohmaver et all . 120071 : iBarret et al.l . 120071 ). Mod- 



els of HFQPOs became more or less popular each 
time that a correlation or no t with the spin was 
found (see Mendez fc Belloriil . 2007 . for a review 
apropos). The original structure of the resonance 
model is unrelated to the spin of the cent ral object . 



It wa s however subsequently proposed l)Lee et al 



20041 ) that in neutron star sources the resonance 



could be excited by direct forcing of the rotating 
star, leading to frequency differences equal to the 
spin or half of it. It is worth noticing that this possi- 
bility is not structural in the resonance model, but 
is a variation of it. On the other hand this brings 
back the problem of the excitation of the modes, 
like in the case of black hole systems. 
Eventually another aspects remains to be explained, 
that is the behavior of the quality factor of HFQPOs 
{Q = J^/FWHM), that accoun ts for the coherence of 
the oscillations. ,Barret et al.l (|2006) reported that 
the quality factor of the lower HFQPOs is higher 
than that of the upper. Furthermore, while the Q of 
the lower frequencies increases with the frequency 
up to a break frequency and then drops, the Q of 
the upper frequencies increases steadily (see Figure 
[3|). These measurements suggest that in different 
sources there is a common physical mechanism not 
only to excite the oscillations, but also to damp 
them. Hopefully this conduct of the quality factor 
may help the theory to to shed some light on these 
mechanism, that for the moment remain obscure. 
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Fig. 3. Measured quality factor of the lower (dots) and upper 
(squares) HFQPOs of 4U 1636-536. The solid and dashed 
lines show the estimates done with the toy model proposed 
in Barret et al. ( 2006i ). from which this figure is taken. 

4. Conclusions 

The nonlinear resonance model for high frequency 
QPOs encounters many internal difhculties. This 
may mark the end of the model or the beginning of 
a new era. 

In our opinion the presence of twin peaks in (al- 
most) commensurate ratio and the detection of sub- 
harmonics still strongly support the idea that weak 
nonlinear resonance is a fundamental ingredient in 
the physics of HFQPOs. 

It is like smelling chocolate in a birthday cake: the 
smell is the unmistakable signature of the presence 
of chocolate. However the other (mathematical and 
physical) ingredients, their relative importance and 
a draft of the recipe (mechanisms) are needed to 
know which cake it is. 
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